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1.  Introduction 


The  past  several  decades  have  seen  growing  interest  in  developing  of  smart  munitions  resulting 
from  attempts  to  increase  accuracy  and  decrease  collateral  damage.  Smart  projectiles  differ  from 
other  guided  weapons,  such  as  guided  missiles,  because  their  electronics  and  control  mechanisms 
must  be  able  to  withstand  extreme  acceleration  loads  associated  with  launch  and  high  spin  rates. 
Control  components  must  be  relatively  inexpensive  since  projectiles  are  typically  fired  in  large 
quantities.  Several  different  types  of  control  mechanisms  have  been  designed  to  meet  these 
requirements — aerodynamic,  thrust,  and  inertial  load  mechanisms.  Common  examples  of 
aerodynamic  mechanisms  are  canards,  gimbaled  nose  configurations,  and  deflection  of  ram  air 
through  side  ports.  Examples  of  thrust  mechanisms  include  cold  gas  jets  and  explosive  thrusters. 
Examples  of  inertial  load  mechanisms  are  rotation  of  an  unbalanced  internal  part  and,  of  specific 
interest  here,  movement  of  an  internal  translating  mass  (ITM).  Oscillation  of  an  internal  mass  at 
the  projectile  roll  frequency  has  been  shown  to  produce  useful  control  authority.  A  design  that 
minimizes  moving  parts  and  substantially  reduces  power  required  is  critical  to  physically 
implement  an  ITM  control  mechanism  on  board  a  smart  munition. 

Previous  investigation  of  projectiles  equipped  with  loose  or  moving  internal  parts  has  revealed 
that  these  configurations  can  result  in  flight  instabilities.  Soper  ( 1 )  considered  the  stability  of  a 
projectile  with  a  cylindrical  mass  fitted  loosely  within  a  cavity.  Using  a  similar  configuration, 
Murphy  (2)  derived  a  quasi-linear  solution  for  the  motion  of  a  projectile  equipped  with  a  moving 
internal  part.  A  detailed  set  of  experiments  was  later  conducted  by  D’Amico  (3)  to  model  the 
motion  of  internal  masses  within  spinning  projectiles  using  a  freely  gimbaled  gyroscope. 

Hodapp  ( 4 )  further  considered  the  effect  of  a  small  offset  between  the  projectile  body  mass  and 
ITM  mass  center.  Hodapp ’s  analysis  of  the  dynamic  equations  for  this  system  showed  that  for 
small  mass  center  offsets,  slight  movement  could  actually  reduce  the  instability  caused  by  the 
loose  internal  part.  Using  controlled  movement  of  an  ITM  as  a  maneuver  control  mechanism 
was  considered  by  Petsopoulous  et  al.  (5)  for  use  on  re-entry  vehicles,  while  Robinett  et  al.  ( 6) 
considered  ITM  control  for  ballistic  rockets.  More  recently,  Menon  et  al.  (7)  examined  ITM 
control  for  use  on  endo-  and  exoatmospheric  interceptors  using  three  orthogonal  ITMs.  Frost 
and  Costello  (3,  9)  have  studied  the  ability  of  an  internal  rotating  mass  unbalance  to  actively 
control  both  fin-  and  spin-stabilized  projectiles.  Most  recently,  Rogers  and  Costello  (10) 
investigated  control  authority  of  a  projectile  equipped  with  a  single  internal  translating  mass.  It 
was  determined  that  significant  control  authority  could  be  created  by  oscillating  the  ITM  at  the 
projectile  roll  frequency. 
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This  report  outlines  a  notional  design  of  an  ITM  actuator  that  generates  sufficient  control 
authority  using  relatively  low  power.  Control  moment  is  generated  due  to  an  axial  drag  offset 
from  the  system  mass  center  caused  by  the  lateral  motion  of  the  ITM.  The  report  begins  with  a 
description  of  the  cantilever  beam  system  that  serves  as  the  translating  mass  and  derives  the 
7-degrees-of- freedom  (DOF)  flight  dynamic  model  used  for  trajectory  predictions.  A  description 
of  the  control  law  and  ITM  electromagnetic  actuators  is  also  provided.  The  dynamic  model  is 
subsequently  employed  to  demonstrate  that  the  cantilever  beam  configuration  provides  sufficient 
control  authority  at  reasonable  power  levels.  Trade  studies  examine  the  effect  of  cantilever  beam 
length  on  force  required  and  battery  size.  The  optimum  natural  frequency  and  damping  ratio  of 
the  beam  are  determined  to  minimize  actuator  control  effort.  Finally,  roll-rate  feedback  control 
is  implemented  to  actively  alter  beam  characteristics  as  the  projectile  roll  rate  changes  during 
flight,  further  decreasing  control  effort  required.  A  final  example  case  demonstrates  that 
sufficient  control  authority  can  be  generated  with  relatively  small  battery  sizes  using  the 
optimum  spring  and  damper  coefficients  and  the  roll-rate  feedback  system. 


2.  Cantilever  Beam  Projectile  Dynamic  Model 


The  cantilever  beam  is  a  fixed-free  elastic  beam,  with  one  end  attached  to  the  projectile  at  point 
L  and  the  free  end  floating  within  the  cavity,  constrained  to  vibrate  in  the  ls  -  J s  plane.  The 
beam’s  first  vibrational  mode  is  the  only  mode  considered  significant  to  dynamic  interaction 
with  the  projectile.  For  this  reason,  the  cantilever  beam  can  be  accurately  modeled  as  a  rigid, 
massless  beam,  with  a  spherical  mass  (considered  to  be  a  pennanent  magnet)  attached  to  the  end. 
A  torsional  spring  and  damper  are  attached  to  the  hinge  point  to  simulate  the  elastic  properties  of 
the  cantilever  beam.  This  dynamically  equivalent  system  is  referred  to  as  the  ITM-Beam.  A 
sketch  of  this  configuration  is  shown  in  figure  1.  Note  that  y  is  defined  as  the  angle  between  the 
ITM-Beam  and  centerline  of  the  projectile.  A  permanent  magnet  is  attached  to  the  end  of  the 
beam  and  can  swing  freely  about  the  hinge.  Force  is  exerted  on  the  magnet  by  electromagnets  on 
both  sides  of  the  cavity  to  move  the  beam  to  a  desired  angle. 
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Figure  1.  The  ITM-Beam  projectile. 


2 


Five  reference  frames  are  used  in  the  development  of  the  equations  of  motion  for  this  system, 
namely  the  inertial,  projectile,  translating  mass,  nonrolling,  and  projectile-fixed  S  reference 
frames.  The  projectile  frame  is  obtained  using  the  standard  aerospace  Euler  angle  sequence  of 
rotations  and  is  related  to  the  inertial  frame  by 
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The  N  frame  is  the  standard  nonrolling  reference  frame  often  used  in  projectile  flight  dynamics 
and  is  defined  by  a  rotation  o  1' -(f)  along  the  I B  axis.  The  S  frame  is  also  fixed  to  the  projectile, 
with  its  origin  at  the  hinge  point.  It  is  defined  so  that  the  ITM-Beam  oscillates  about  the  Ks  axis 
and  Js  points  to  the  rear  of  the  cavity  exactly  equidistant  from  both  electromagnetic  actuators. 
Therefore,  the  S  frame  can  be  related  to  the  B  frame  by  two  constant  Euler  angles,  y/s  and  0S ,  so 
that 
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Throughout  the  rest  of  this  report,  fixed  angles  of  y/s  =  90°  and  Bs  =  0  are  used;  thus,  the  S  frame 
can  be  obtained  by  a  single  90°  rotation  about  the  KB  axis,  resulting  in  the  orientation  shown  in 
figure  1 .  The  T  frame  is  fixed  to  the  ITM-Beam  and  is  related  to  the  S  frame  by  the  relationship 
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Note  that  the  T  frame  is  aligned  with  the  S  frame  when  y  =  0  .  All  equations  in  this  report  use  the 
following  shorthand  notation  for  trigonometric  sine,  cosine,  and  tangent  functions:  sa  =  sin  a , 
ca  =  cos  a ,  and  ta  =  tan  a  . 

Throughout  the  development  of  the  equations  of  motion,  two  operators  will  be  used:  C  to 
denote  components  of  a  vector  in  a  specific  frame  and  S  to  denote  the  skew  symmetric  cross 
product  operator.  The  vector  component  operator  outputs  a  column  vector  comprised  of  the 
components  of  an  input  vector  in  a  given  frame.  For  example,  if  the  position  vector  from  a  to  ft 
is  expressed  in  reference  frame  A  as  ra^[s  =  Axa/jI  a  +  AyafjJ a  +  Azu/jK  a  ,  then  the  vector 
component  operator  acting  on  this  vector  yields 
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C  A(ra^p)  =  \ 


to, aP 

A yaP 

to. aP 


(4) 


Notice  that  the  reference  frame  is  denoted  by  the  subscript  on  the  operator. 


The  cross-product  operator  outputs  a  skew  symmetric  matrix  using  the  components  of  an  input 
vector  in  the  reference  frame  denoted  in  the  subscript.  For  example,  if  the  position  vector  from  a 
to  f>  is  expressed  in  reference  framed  as  ra^fj  =  AxufjI  a  +  AvafjJ  a  +  Azu/jK  a  ,  then  the  cross-product 
operator  acting  on  ra^f,  expressed  in  reference  frame  A  is 
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2,1  Kinematics 

The  velocity  of  the  composite  body  mass  center  can  be  described  in  the  inertial  frame  or  the 
projectile  reference  frame. 

vcu  =  xl i  +  yJ  i  +  zKi  =  uIb  +  vJ  b  +  wK  «  .  (6) 


The  translational  kinematic  differential  equations  relate  these  two  representations  of  the  mass 
center  velocity  components. 
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The  angular  velocity  of  the  projectile  with  respect  to  the  inertial  reference  frame  can  be  written 
in  terms  of  appropriate  Euler  angle  time  derivatives  or  in  terms  of  projectile  frame  angular 
velocity  components. 

<X>BH  =  <j)I  B  +  OJ  N  +  \jfKl  —  pi  B  +  CjJ  B  +  f  Kb  .  (8) 


The  kinematic  relationship  between  time  derivatives  of  the  Euler  angles  and  projectile  reference 
frame  angular  velocity  components  represents  the  rotational  kinematic  differential  equations. 
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The  final  kinematic  differential  equation  is  the  trivial  relationship 

r=^Beam-  (N>) 

2.2  Dynamics 

The  translational  dynamic  equations  for  the  ITM-Beam  projectile  are  derived  through  force 
balancing.  Force-balance  equations  for  the  projectile  and  ITM-Beam  are  given  as 

mPaPII  =Wp+Fp-Fc-Fj  (11) 

and 

ntTaXII  =WT  +  Fc  +  Fj  (12) 

respectively,  where  mp  and  mT  are  the  masses  of  the  projectile  and  ITM-Beam,  ap/l  and  axn  are 
the  accelerations  of  points  P  and  X  with  respect  to  the  inertial  frame,  and  Wp,  WT ,  FP,FI  and  F 
are  the  weights  of  the  projectile  and  ITM-Beam,  total  aerodynamic  force  exerted  on  the 
projectile,  input  force  exerted  on  the  ITM-Beam  by  the  actuators,  and  hinge  constraint  force, 
respectively.  The  definition  of  the  system  center  of  mass  leads  to 

maCII  =mpaPll  +mTax/1  .  (13) 

Therefore,  adding  equations  1 1  and  12  and  noting  the  relationship  in  equation  13,  the 
translational  dynamic  equation  for  the  system  is  formed. 

macn  =  WP  +  WT  +  Fp  .  (14) 

The  aerodynamic  forces  given  by  Fp  in  equation  14  are  obtained  using  the  standard  aerodynamic 
expansion  employed  for  projectile  flight  dynamic  simulation.  Both  steady  aerodynamic  forces 
and  Magnus  forces  are  included  as  well  as  steady  and  unsteady  aerodynamic  moments.  The 
aerodynamic  coefficients  and  aerodynamic  center  distances  used  to  generate  these  forces  and 
moments  are  all  a  function  of  the  local  Mach  number  at  the  center  of  mass  of  the  projectile. 
Computationally,  these  Mach  number-dependent  parameters  are  obtained  by  a  table  look-up 
scheme  using  linear  interpolation.  A  full  description  of  the  weight  force  and  body  aerodynamic 
forces  and  moments  is  provided  in  reference  10.  Writing  equation  14  in  the  projectile  reference 
frame  yields 
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Note  that  XB,YB,  and  ZB  are  projectile  reference  frame  components  of  the  sum  of  the  three 
forces  given  in  equation  14. 


The  rotational  dynamic  equations  are  obtained  by  first  equating  the  /  frame  time  rate  of  change 
of  the  system  angular  momentum  about  the  system  mass  center  to  the  total  applied  external 
moments  on  the  system  about  the  system  mass  center  in  the  Is  and  Js  directions,  given  by 
equations  16  and  17.  Then,  the  same  moment  equation  is  used  for  each  body  separately,  this 
time  written  in  the  direction.  These  four  equations  are  given  by 
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(17) 

1{WT}  +  kTy  +  kDy  , 

(18) 

and 


r  i 


Kc 


dH 


B/I 


dt 


+  r, 


L^P 


\ 

x  mpdpn 

) 


=  f„p„BLcr  +KS.  S(Fw,)Cs  (WF)-kTr-kDf  +  Ks 

•CS(I Mi), 


(19) 


where,  as  before,  cy  =  cos  (7) .  Note  that  equation  18  equates  the  time  rate  of  change  of  the 
angular  momentum  of  the  ITM-Beam  to  the  total  moment  applied  to  the  ITM-Beam.  Equation 
19  equates  the  time  rate  of  change  of  the  angular  momentum  of  the  projectile  to  the  total  moment 
applied  to  the  projectile. 

Several  intermediate  expressions  will  be  useful  in  deriving  the  rotational  dynamic  equations  in 
the  body-fixed  S  frame.  First,  note  that  the  well-known  two  points  fixed  on  a  rigid  body  formula 
yields  the  relationship 


^xii  UCII  +  [^b/j  xi^i +<5r/J  xi^  +ryB/J  x(ryfi// xfp  )  +  (yr// x(hir/J  xi^  )J  .  (20) 

m  L  J 


Equation  20  is  used  to  expand  dx/I  in  terms  of  known  quantities  and  state  derivatives.  Also, 
using  the  definition  of  the  system  center  of  mass,  it  can  be  shown  through  algebraic  manipulation 
that 


rL^>P  X  mpaptl  +  rL^X  X  mTaXII  rL^>P  X  maC/I  +  fp->AT  X  mTClx/I  . 


(21) 
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Equation  21  is  also  used  to  expand  the  cross-product  terms  on  the  left-hand  side  of  equations  16 
and  17  in  terms  of  known  quantities  and  state  derivatives. 

Equations  16-19  can  be  expanded  using  the  expressions  in  equations  20  and  21  and  rearranged 
to  fonn  a  4  x  4  system  of  equations  given  by 


"4i 

Ai 

As 

A\  4 

p 

A2i 

A22 

A.s 

A4 

q 

A 

A, 

A32 

As 

Aa 

r 

Bs 

_A4\ 

A42 

As 

- 1 

”3" 

j\ 

B4_ 

In  equation  22,  rows  1-4  correspond  to  equations  18,  19,  17,  and  16,  respectively.  The  full 
expressions  for  the  values  of  the  A  matrix  and  B  vector  are  lengthy  and  provided  in  appendix  A. 
The  set  of  equations  given  by  7,  9,  10,  15,  and  22  constitutes  the  equations  of  motion  for  the  ITM- 
Beam  projectile.  Given  a  known  set  of  initial  conditions,  these  14  scalar  equations  are 
numerically  integrated  forward  in  time  using  a  fourth-order  Runge-Kutta  algorithm  to  obtain  a 
single  trajectory. 


2.3  Description  of  Controller 

To  create  trajectory  alterations,  the  ITM-Beam  must  be  moved  in  a  prescribed  manner.  The 
control  law  is  formulated  based  on  a  feedback  linearization  technique  (11),  which  assumes  full- 
state  feedback.  To  compute  the  control  force,  finut ,  required  to  deflect  the  ITM-Beam  to  the 
desired  angle,  use  the  following: 


f.  = 

J  input 


-1 

BLc 


^  ilc  AuP  AyA  A\A‘  AA  j  Kjt 


command 


~K  (  y  (y  _  y) 

1  \/  command  /  /  2  \/  command  /  /’ 


(23) 


Note  that  Bflc  is  defined  in  appendix  A  and  derived  from  equation  18.  Likewise, 

An,  Ai2 ,  Al3 ,  and  An  are  from  equation  1 8  and  provided  in  appendix  A.  Note  that  equation  1 8  is 
used  to  compute  the  feedback  linearization  control  rather  than  equation  19  since  equation  19 
would  require  feedback  of  aerodynamic  loads.  This  is  a  complicated  task  and  can  be  avoided  by 
using  equation  18  instead. 


The  commanded  deflection  angle,  ycommani ,  is  generated  by  synchronizing  ITM-Beam  movement 
with  the  projectile  roll  angle.  This  is  done  by  setting 


(A 

T command  =  Sil1  '  C°s(0+&) 

V  IJLs 


(24) 


where  A  is  the  magnitude  of  oscillation  of  point  X  from  the  cavity  center  and  <j>T  is  a  trim  angle 
used  to  define  the  plane  of  control.  Derivatives  of  equation  24  are  computed  analytically  and 
used  in  equation  23.  Note  that  the  feedback  linearization  controller  is  developed  and  used 
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within  the  simulation  solely  to  create  the  prescribed  motion  of  the  ITM-Beam  for  control 
authority  analysis  and  determination  of  power  requirements.  In  this  case,  its  purpose  is  to  match 
the  ITM-Beam  oscillation  frequency  to  the  projectile  roll  rate  within  the  simulation. 

2.4  Description  of  Electromagnetic  Actuator  Control  System 

A  zoom  view  of  the  ITM-Beam  mechanism  is  shown  in  figure  2.  Two  electromagnets,  each  at 
opposite  ends  of  the  cavity,  exert  force  on  the  fixed  magnet  at  the  end  of  the  ITM-Beam. 


Figure  2.  Zoom  view  of  the  ITM-Beam  system. 


The  force  exerted  on  a  fixed  magnetic  dipole  DM  is  given  by  Purcell  (12)  as 

dB 

finpu,=DM — (25) 

The  magnetic  dipole  moment  per  unit  mass  is  a  unique  property  of  a  material,  with  units 
joules/tesla/slug.  For  example,  assuming  the  ITM  is  made  of  magnetized  iron  and  using  an  ITM- 
Beam  mass  of  mt  =  0.05  slugs,  the  dipole  moment  is  found  to  be  DM=  171.5  J/T.  The  quantity 
dB 

— -  can  be  found  by  first  recognizing  that  the  magnetic  field,  B  ,  of  an  iron  core  solenoid  is 
given  by  the  expression 

B  _ 
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where  IE  is  the  current  through  the  electromagnet,  n  is  the  magnetic  constant  (4ttx  10  N/A2),  zM 
is  the  distance  from  the  endpoint  of  the  ITM  beam  to  the  nearest  electromagnet  actuator,  k  is  the 
dimensionless  relative  permeability  of  iron  (200  at  a  magnetic  flux  density  of  0.002  W/m  ),  n  is 
the  number  of  coils  per  meter,  b  is  the  radius  of  the  solenoid,  and  La  is  the  length  of  the  solenoid. 
For  all  the  following  cases,  values  of  3,  2,  and  10000  cm  were  used  for  b,  L,  and  n,  respectively. 
Taking  the  derivative  of  equation  26  with  respect  to  zM , 


dB  IE/ukn 

d*M  2 


(27) 


At  each  timestep,  control  force  is  computed  using  the  feedback  linearization.  Knowing  the 

dB 

magnetic  dipole  moment  and  required  control  force,  the  quantity  — -  is  computed  at  each 

SzM 

timestep  using  equation  25.  Then,  knowing  the  position  of  the  ITM-Beam  with  respect  to  the 
actuators,  the  current  required  can  be  computed  at  each  timestep  by  rearranging  equation  27  so 
that 


/ 


E 


1-1 


SB 

_ z_ 


(28) 


Note  that  the  electromagnet  dimensions  used  in  equation  28  are  comparable  to  dimensions  for 
commercially  available  iron-core  electromagnets. 

2.5  Description  of  Example  Projectile 

The  example  projectile  used  in  the  following  simulations  is  a  representative  fin-stabilized 
projectile.  Relevant  example  projectile  and  ITM-Beam  dimensional  and  mass  properties  are 
outlined  in  table  1.  The  hinge  point  is  0.9  ft  behind  the  projectile  mass  center  ( P ),  and  the  ITM- 
Beam  oscillation  amplitude  is  given  by  SLsin(ymax)  =  0.157  ft,  unless  otherwise  specified.  In 

the  following  cases,  the  projectile  is  traveling  through  a  standard  atmosphere,  with  no 
atmospheric  wind. 

Table  1.  Relevant  example  projectile  and  ITM-Beam  properties. 


Projectile  Mass  (slugs) 

1.206 

ITM-Beam  Mass  (slugs) 

0.05 

Projectile  Reference  Diameter  (ft) 

0.344 

Projectile  Mass  Center  Position 
Measured  Along  Stationline  (ft) 

1.18 

Projectile  Roll  Inertia  (slugs-ft2) 

0.0278 

Projectile  Pitch  Inertia  (slugs-ft2) 

0.6291 

ITM-Beam  Length  (ft)  (BL) 

0.3 
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3.  Results 


An  example  trajectory  of  the  ITM-Beam  projectile  is  compared  to  an  example  trajectory  of  a 
projectile  equipped  with  a  strictly  translating  internal  mass  for  model  validation  purposes.  The 
translating  mass  projectile’s  dynamic  equations  are  given  in  reference  10,  and  a  previously 
validated  model  of  this  system  was  used  for  trajectory  predictions.  A  diagram  of  the  translating 
mass  projectile  is  provided  in  appendix  B.  Initial  conditions  used  for  the  example  trajectory 
were  x  =  0.0  ft,  y  =  0.0  ft,  z  =  0.0  ft,  u  =  282 1 .0  ft/s,  v  =  0.0  ft/s,  w  =  0.0  ft/s,  <f>  =  1 .5707  rad,  6 
=  0.05  rad,  y/  =  0.0  rad ,p  =  5.0  rad/s,  q  =  0.0  rad/s,  and  r  =  0.0  rad/s.  The  feedback  linearization 
gains  were/f0  =  0.001,  Kx  =  1000.0,  and  K2  =  5  x  105.  In  the  following  cases,  the  ITM-Beam 
oscillation  frequency  is  locked  to  the  projectile  roll  rate.  Figures  3-6  show  the  trajectories  for 
the  ITM-Beam  projectile,  translating  mass  projectile,  and  the  rigid  projectile  with  no  internal 
moving  mass  (denoted  “Rigid  6DOF”).  The  two  translating  mass  trajectories  are  generated 
solely  to  demonstrate  control  authority  and  validate  the  ITM-Beam  simulation;  thus,  both 
controlled  rounds  are  commanded  to  maximum  possible  deflection.  Notice  that  the  trajectories 
of  the  internal  mass  projectiles  match  nearly  identically,  even  though  the  dynamic  equations  for 
the  two  systems  are  significantly  different.  The  correlation  between  these  two  models  serves  as 
validation  of  the  ITM-Beam  simulation. 


Figure  3.  Altitude  vs.  range  for  example  trajectory. 
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Figure  5.  Velocity  (u)  vs.  time  for  example  trajectory. 
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Time  (s) 


Figure  6.  Roll  rate  vs.  time  for  example  trajectory.  The  thick  lines  represent  high- 

frequency  oscillations,  which  occur  only  for  the  ITM-Beam  and  translating 
mass  case. 


Figure  7  shows  a  selected  time  history  of  ITM  displacement  from  the  projectile  centerline.  For 
the  ITM-Beam  projectile,  this  displacement,  sx ,  is  given  by 

sx=BLsin(y)  .  (29) 

Notice  that  the  two  time  histories  for  ITM  displacement  shown  in  figure  7  are  nearly  identical.  A 
control-force  time  history  for  an  example  ITM-Beam  projectile  simulation  can  be  used  together 
with  an  ITM  displacement  time  history  to  generate  a  time  history  of  current  required  for  each 
electromagnetic  actuator.  This  is  accomplished  using  the  procedure  outlined  in  equations  25-28. 
Furthermore,  this  current  time  history  can  be  integrated  to  produce  the  total  charge  required  for  a 
given  example  flight  in  A/s.  This  value  for  total  charge  can  be  used  to  size  the  battery  for  the 
ITM-Beam  control  system.  The  total  charge  required  for  the  example  flight  shown  previously 
was  13.5  A/s.  Note  that  this  number  is  relatively  large  since  the  spring  and  damper  coefficient 
values  have  not  been  optimized  for  this  preliminary  example  case. 
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Figure  7.  Selected  time  histoiy  of  ITM  displacement  from  projectile  centerline. 


Figure  8  shows  a  segment  of  the  current  time  history  for  the  example  ITM-Beam  simulation  used 
previously.  Notice  that  f  (labeled  “Force”  in  the  plot)  and  the  ITM  displacement,  sx ,  are 
shown  on  the  same  plot  as  the  current  time  history  to  demonstrate  the  phase  relationships 
between  current,  input  force,  and  ITM  displacement.  When  the  ITM-Beam  is  displaced  in  the 
positive  /  direction  (sY  positive),  electromagnet  EM  1  shuts  off,  and  EM  2  is  responsible  for 
control.  Likewise,  when  the  ITM-Beam  is  displaced  in  the  negative  /  direction  ( sx  negative), 
electromagnet  EM  2  shuts  off,  and  EM  1  is  responsible  for  control.  This  scheme  takes  advantage 
of  the  fact  that  the  electromagnets  are  more  effective  when  the  ITM  is  at  close  range.  Since 
current  required  is  a  nonlinear  function  of  distance  to  the  ITM  and  control  force  required,  the 
current  time  history  is  not  sinusoidal  like  the  ITM  displacement  and  control-force  time  histories. 

The  length  of  the  ITM-Beam  has  a  significant  impact  on  the  force  required  to  move  the  beam  in 
a  prescribed  fashion.  From  equation  18,  the  external  moment  exerted  on  the  ITM-Beam  by  the 
actuators  in  the  K s  direction  about  point  L  is  given  by 

Ks-MLBeam=~flnputBLcr  .  (30) 
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Figure  8.  Segment  of  current  vs.  time  for  ITM-Beam  actuators. 


Therefore,  the  input  force  required  to  exert  a  given  control  moment  on  the  ITM-Beam  varies 
inversely  with  beam  length.  In  addition,  for  a  given  maximum  ITM  displacement,  sx  ,  the 
maximum  angular  displacement,  ymax ,  varies  inversely  with  beam  length.  Trade  studies  verified 
these  results  using  example  simulations  with  various  beam  lengths.  The  maximum  ITM 
displacement  from  the  projectile  centerline  was  sx  =  0.157  ft  for  all  cases.  Figures  9  and  10 
show  maximum  angular  displacement  and  average  force  required,  respectively,  as  a  function  of 
beam  length.  Notice  that  as  beam  length  increases,  maximum  angular  displacement  and  average 
force  required  decrease.  The  average  control-force  reduction  with  increased  beam  length  occurs 
due  to  the  increased  “efficiency”  of  the  ITM-Beam  actuators  (i.e.,  the  same  control  moment 
requires  less  control  force  for  a  larger  beam  length). 
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Figure  9.  Maximum  angular  displacement  vs.  beam  length. 


Figure  10.  Average  force  required  vs.  beam  length. 
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As  outlined  previously,  the  ITM-Beam  system  is  a  dynamic  model  used  to  represent  a  fixed-free 
elastic  beam.  The  rigid  ITM-Beam  is  attached  to  the  projectile  at  point  L  with  a  torsional  spring 
and  torsional  damper  to  model  the  elastic  beam’s  vibrational  properties.  A  trade  study  examines 
how  force  and  power  requirements  vary  with  different  spring  and  damping  coefficients.  Once  an 
optimum  spring  and  damper  coefficient  are  determined,  they  can  be  used  to  identify  the  proper 
elastic  properties  of  the  fixed-free  beam  for  a  prototype  system. 

The  performance  of  the  system  is  examined  for  a  range  of  torsional  spring  constants  and 
damping  ratios  for  the  example  projectile  rolling  at  a  steady-state  rate  of -128  rad/s.  The 
projectile  trajectory  is  simulated  for  a  2-s  flight  with  no  gravity.  This  simplified  flight  profile  is 
used  solely  to  establish  the  correlation  between  spring  and  damper  parameters  and  average  force, 
average  power,  and  total  battery  charge  required.  Figure  1 1  shows  the  projectile  roll-rate  time 
history  for  this  flight  profile.  The  high-frequency  oscillation  of  the  roll  rate  occurs  at  the  mass 
oscillation  frequency.  This  is  due  to  the  continually  changing  axial  moment  of  inertia  of  the 
projectile  as  the  mass  translates.  When  the  translating  mass  travels  farther  from  the  centerline, 
axial  inertia  grows  and  roll  rate  decreases  due  to  conservation  of  angular  momentum.  Likewise, 
when  the  translating  mass  returns  toward  the  centerline,  axial  inertia  decreases  and  the  roll  rate 
increases. 


Figure  11.  Roll  rate  vs.  time  for  example  simulation. 
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Figures  12-14  show  the  effect  of  spring  and  damper  coefficients  on  average  force,  average 
power,  and  total  battery  charge  required.  In  figure  12,  it  can  be  clearly  seen  that  an  optimum 
torsional  spring  constant  exists  in  which  the  spring-mass-damper  system  of  the  ITM-Beam 
operates  near  resonance  with  the  projectile  roll  rate.  These  peaks  are  not  as  sharp  as  typical 
spring-mass-damper  resonant  peaks  due  to  the  fact  that  the  projectile  roll  rate  varies  over  time. 
However,  note  that  significant  reductions  in  force  are  achieved  if  the  spring  constant  is  placed 
near  its  optimum  value  and  damping  is  lowered  as  much  as  possible.  These  reductions  in  force 
are  mirrored  by  reductions  in  average  power  and  total  charge  required,  resulting  in  significantly 
smaller  battery  sizes. 

A  similar  study  examines  the  same  spring-mass-damper  parameters  for  a  full  flight  profile  of  the 
example  projectile  using  the  same  initial  conditions  as  those  used  in  the  first  example  study. 
Figures  15-17  show  that  as  in  the  partial  flight  profile  case,  optimal  spring  coefficients  can  be 
found. 


Figure  12.  Average  force  required  vs.  torsional  spring  constant,  partial  flight  profile. 
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Figure  13.  Average  power  required  vs.  torsional  spring  constant,  partial  flight  profile. 


Figure  14.  Total  charge  required  vs.  torsional  spring  constant,  partial  flight 
profile. 
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Figure  15.  Average  force  required  vs.  torsional  spring  constant,  frill  flight 
profile. 


Figure  16.  Average  power  required  vs.  torsional  spring  constant,  full  flight 
profile. 
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Figure  17.  Total  charge  required  vs.  torsional  spring  constant,  full  flight  profile. 


However,  the  results  for  the  partial  flight  profile  have  a  significantly  sharper  peak  than  the  results 
for  the  full  flight  profile.  This  is  because,  as  shown  in  figure  6,  roll  rate  of  the  projectile  varies 
between  5  rad/s  initially  and  a  final  value  of  ~80  rad/s.  This  large  variation  in  roll  rate  means 
that  the  spring  coefficient  is  only  optimized  for  a  very  short  period  of  the  overall  flight,  and  the 
broad  peaks  shown  in  figures  15-17  result.  To  demonstrate  this,  figure  18  shows  a  current  time 
history  for  an  example  full  flight  trajectory  using  kT  =  70.0  lb/rad  and  kD  =  0.05  lb/rad/s.  Note 
that  the  spring  constant  is  optimized  for  the  projectile  roll  rate  ~1  s  into  flight  and  once  again 
after  spin  decay  occurs  ~7  s  into  the  flight. 

Despite  the  broad  nature  of  the  peaks  shown  in  figures  15-17,  significant  size  and  weight  savings 
can  be  achieved  using  the  proper  spring  constants  in  the  form  of  smaller  batteries.  As  shown  in 
figure  17,  batteries  with  a  total  charge  of  less  than  5  A/s  may  be  used  for  systems  with  optimal 
spring  coefficients  and  low  damping  ratios.  Furthermore,  figure  18  shows  that  reasonable 
maximum  current  levels,  on  the  order  of  150  mA,  can  be  expected  with  an  optimized  system. 
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Figure  18.  Current  through  actuators  vs.  time  for  example  full  flight 
trajectory. 


Average  force  levels  and,  therefore,  total  charge  required  can  be  decreased  even  further  by 
actively  changing  the  elastic  properties  of  the  beam  during  flight.  A  fixed-free  cantilever  beam 
like  that  used  in  this  system  has  a  first  vibrational  mode  shape  of 


sx 


=  sin 
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where  xB  is  the  distance  along  the  beam.  The  natural  frequency  of  the  first  vibrational  mode  is 


K  l~E 
2  BL  ]]  p  ’ 


(32) 


where  E  is  the  modulus  of  elasticity  of  the  material  and  p  is  the  density  of  the  beam  (13).  By 
changing  the  modulus  of  elasticity,  it  is  therefore  possible  to  tune  the  natural  frequency  of  the 
cantilever  beam  to  a  desired  value.  Recent  investigations  into  smart  materials  (14-17), 
specifically,  materials  used  in  tunable  vibration  absorbers,  have  shown  that  various  methods  can 
be  used  to  actively  alter  a  material’s  modulus  of  elasticity,  allowing  the  beam’s  torsional  spring 
constant  to  be  actively  optimized  during  flight  as  the  projectile  roll  rate  changes.  This  would 
allow  the  ITM-Beam  system  to  operate  with  the  lowest  possible  power  through  the  entire  flight, 
yielding  further  reductions  in  battery  size. 
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To  investigate  this,  several  example  simulations  were  run.  The  first  set  simulated  the  projectile 
for  the  full  flight  using  the  optimum  spring  constants  obtained  from  figure  17.  This  produced  the 
least  possible  battery  charge  required  for  the  ITM-Beam  system,  with  a  fixed  spring  constant  for 
each  damping  ratio  considered.  The  second  set  of  simulations  included  a  roll-rate  feedback 
mechanism.  In  these  cases,  at  specific  points  throughout  the  flight,  the  torsional  spring  and 
damping  coefficients  of  the  ITM-Beam  were  adjusted  to  match  the  roll  rate.  Figure  19  shows 
how  the  torsional  spring  constant  was  adjusted  for  an  example  flight  with  C=  0.05.  Note  that  the 
curve  in  figure  19  has  the  same  qualitative  shape  as  the  roll  rate  time  history  shown  in  figure  6. 
Table  2  summarizes  the  results  of  the  two  sets  of  simulations  with  the  fixed  torsional  spring  and 
variable  torsional  spring  constants.  Note  that  in  the  variable  torsional  spring  cases,  the  torsional 
damping  coefficient  kD  was  adjusted  slightly  to  keep  the  damping  ratio  C  constant.  Figure  20 
demonstrates  that  implementing  the  roll-rate  feedback  system  saves  ~1  A/s  of  charge  for  all 
damping  ratios  considered  (34.4%  decrease  in  battery  size  for  C=  0.02,  28.1%  decrease  for  C 
=  0.05,  20.4%  for  C=  0.1,  and  1 1.2%  for  C=  0.2). 


Figure  19.  Torsional  spring  constant  vs.  time  for  roll-rate  feedback  system,  f  =  0.05. 
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A  conceptual  design  for  a  projectile  equipped  with  a  controllable  internal  mass  is  proposed  using 
a  cantilever  beam  configuration.  The  system  is  designed  to  minimize  moving  parts  and  exploit 
the  beam’s  elastic  properties,  thereby  creating  a  robust  mechanism  that  minimized  power 
requirements.  A  dynamic  model  of  the  configuration  is  developed,  as  is  a  model  of  the 
electromagnetic  actuator  system  to  study  current  requirements  and  battery  size.  Example 
simulations  show  that  useful  control  authority  is  generated  with  this  projectile  control 
mechanism,  matching  the  results  from  earlier  studies  of  projectiles  with  internal  translating  mass 
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control.  Furthermore,  by  tuning  the  natural  frequency  of  the  vibrating  beam  to  the  projectile  roll 
rate,  significant  reductions  in  power  requirements  are  achieved  on  the  order  of  60%.  Power 
reductions  can  be  further  realized  with  active  tuning  of  the  beam’s  elastic  properties  during 
flight.  Battery  size  requirements  to  power  the  actuator  are  modest,  with  several  commercial-off- 
the-shelf  options  available.  The  internal  oscillating  beam  configuration  shows  promise  as  a 
viable,  cost-effective,  reliable  projectile  control  mechanism. 
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The  terms  of  equation  22  are  given  as  follows: 
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Figure  B-l.  Translating  mass  projectile  schematic. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


flj ,  a2 ,  «3  Components  in  the  S  frame  of  the  acceleration  of  the  system  mass  center  with 

respect  to  the  inertial  frame 

aPII  Acceleration  of  the  projectile  mass  center  with  respect  to  the  inertial  frame 

a X/I  Acceleration  of  the  end  of  the  ITM-Beam  with  respect  to  the  inertial  frame 

A  Maximum  magnitude  of  ITM-Beam  displacement  from  center  of  cavity 

b  Radius  of  the  electromagnetic  actuator 

BL  Length  of  the  ITM-Beam  assembly 

B_  Magnetic  field  produced  by  the  electromagnetic  actuators  in  the  I s  direction 

C  Composite  body  center  of  mass 

DM  Magnitude  of  the  magnetic  dipole  moment  of  the  fixed  magnet  at  the  end  of  the 

ITM-Beam 

Fc  Hinge  constraint  force  on  the  internal  translating  mass 

F1  Input  force  exerted  by  electromagnetic  actuators  on  ITM-Beam 

Fp  Total  aerodynamic  force  exerted  on  the  projectile 

fi  t  Scalar  value  of  the  total  force  applied  to  ITM-Beam  by  electromagnetic  actuators 

HB/I  Angular  momentum  of  the  projectile  with  respect  to  the  inertial  frame  about  the 

projectile  center  of  mass 

HPII  Angular  momentum  of  the  ITM-Beam  with  respect  to  the  inertial  frame  about 

point  X 

fTxx  The  xx  component  of  the  ITM-Beam  moment  of  inertia  matrix  about  point  X 

IPxx  The  xx  component  of  the  projectile  moment  of  inertia  matrix 

KI  Inertial  frame  unit  vectors 
I B,JB,  Kb  Projectile  reference  frame  unit  vectors 


IN,  JN,  Kn  Nonrolling  reference  frame  unit  vectors 
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1 "y  ^  J~ p  ^  p 


Z 

44 

ntp 

mT 

m 

Mi 


ml 

system 

mL_ 


'  L^P 


ITM-Beam  fixed  frame  unit  vectors 
S  frame  unit  vectors 

Current  through  the  electromagnetic  actuator 
Torsional  damper  coefficient 
Torsional  spring  coefficient 

Junction  of  between  the  ITM-Beam  and  the  projectile,  referred  to  as  the  “hinge 
point” 

Length  of  the  electromagnetic  actuator 
Mass  of  the  projectile  with  the  cavity 

Mass  of  the  ITM-Beam 
Total  mass  of  the  system 

External  moments  applied  to  the  projectile  about  the  hinge  point  L 

External  moments  applied  to  the  projectile-ITM  system  about  the  hinge  point  L 

External  moment  exerted  on  the  ITM-Beam  by  the  actuators  in  the  Ks  direction 
about  the  hinge  point  L 

Projectile  center  of  mass 

Distance  vector  from  the  hinge  point  L  to  point  X  at  the  end  of  the  ITM-Beam 

Distance  vector  from  the  hinge  point  L  to  the  projectile  center  of  mass  P 

Distance  vector  from  projectile  center  of  mass  P  to  point  X  at  the  end  of  the  ITM- 
Beam 


Components  in  the  S  frame  of 


V 

r4,  r5,  r6  Components  in  the  S  frame  of  p^x 

p,  q,  f  Components  of  cbBII  in  the  S  frame 

p ,  q,  r  Components  of  coBII  in  the  projectile  reference  frame 
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U,  V,  w 


WP 

WT 

X 

X,  y,  z 


Y  Y  7 

^  B>  1B>  ^ B 


ccTu 

Y 

0,  y/ 


Oj,  y/ P 


co, 


BII 


(Or, 


T!I 


CO, 


Beam 


Translational  velocity  components  of  the  composite  body  center  of  mass  resolved 
in  the  projectile  reference  frame 

Velocity  of  the  system  mass  center  with  respect  to  the  inertial  frame 
Weight  of  the  projectile  (without  the  ITM-Beam) 

Weight  of  the  ITM-Beam 
Point  at  the  end  of  the  ITM-Beam 

Position  vector  components  of  the  composite  body  center  of  mass  expressed  in 
the  inertial  reference  frame 

Total  external  force  components  on  the  projectile  and  ITM-Beam  system 
expressed  in  the  projectile  reference  frame 

Angular  acceleration  of  the  projectile  body  with  respect  to  the  inertial  frame 
Angular  acceleration  of  the  ITM-Beam  with  respect  to  the  inertial  frame 
Deflection  angle  of  the  ITM-Beam 
Euler  roll,  pitch,  and  yaw  angles 

Euler  pitch  and  yaw  angles  for  the  orientation  of  the  S  frame  with  respect  to  the  B 
frame 

Angular  velocity  of  the  projectile  body  with  respect  to  the  inertial  frame 
Angular  velocity  of  the  ITM-Beam  with  respect  to  the  inertial  frame 
Magnitude  of  the  angular  velocity  of  the  ITM-Beam  with  respect  to  the  S  frame 
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